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ABSTRACT 

We present a new observational method to type the explosions of young supernova remnants (SNRs). 
By measuring the morphology of the Chandra X-ray line emission in seventeen Galactic and Large 
Magellanic Cloud SNRs with a multipole expansion analysis (using power ratios), we find that the 
core-collapse SNRs are statistically more asymmetric than the Type la SNRs. We show that the two 
classes of supernovae can be separated naturally using this technique because X-ray line morphologies 
reflect the distinct explosion mechanisms and structure of the circumstellar material. These findings 
are consistent with recent spectropolarimetry results showing that core-collapse SNe are intrinsically 
more asymmetric. 

Subject headings: methods: data analysis — supernovae: general — supernova remnants — X-rays: 
ISM 



X 

O : 

to: 



> 

00 

o 

m 



On 
O 



1. INTRODUCTION 

Over the past few decades, observations and theory of 
supernovae explosions (SNe) and their remnants have ad- 
vanced tremendously. The original classification of SNe 
was based on their optical spectra near maximum bright- 
ness (Minkowski 1941; see Filippenko 1997 for a review), 
when features originate from the outermost layers of the 
exploding stars. At later times, SN expansion causes the 
chemical and physical structure of the SN ejecta and the 
circumstellar material (CSM) to dominate the spectra. 
Study of these nebular spectra revealed two physically 
disinct classes of SNe: Type la and core-collapse (CC; 
Types lb, Ic, and II) explosions. There is a broad consen- 
sus now that Type la SNe constitute the thermonuclear 
detonation of a white dwarf, while CC SNe follow the 
gravitational collapse of a massive (> 8M ) star. These 
different pathways produce distinct observational signa- 
tures that reflect their underlying explosion mechanisms. 

As the ejecta continue to expand following an explo- 
sion, they are heated to X-ray emitting temperatures 
by collisionless shocks. Due to this phenomenon, X-ray 
emission lines are observable in young supernova rem- 
nants (SNRs with ages < 10,000 years) from the metals 
produced within the stellar interiors and during the SNe. 
Thus, study of X-ray spectra and X-ray morphology can 
reveal much about SNe progenitors (see Weisskopf & 
Hughes 2006 for a review). For example, the strength 
of X-ray emission lines can be used to infer chemical 
abundances and to compare with the values predicted 
by Type la and CC SNe models. In particular, the O/Fe 
ratio is often used for this purpose, since Type la SNe 
produce comparatively more Fe than CC SNe, and CC 
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SNe yield large amounts of O. However, the analysis of 
eject a-dominated X-ray spectra is very challenging given 
our limited understanding of a number of topics, includ- 
ing the hydrodynamic interaction between the ejecta and 
the ambient medium, the inner structure of the ejecta, 
and the electron heating processes at the reverse shock 
(Badenes et al. 2003, Rakowski et al. 2006). So far, 
this detailed analysis of the ejecta emission including all 
these subtleties has only been done for a few objects. 

In this Letter, we present a new, observational method 
for the typing of young SNRs. Specifically, we describe 
how the asymmetry of X-ray line emission from SNRs 
can be used to discern between a Type la or a CC ori- 
gin. Our technique is general and quantitative, and it 
is independent of the conditions of the plasma. Thus, 
this approach can be used systematically for comparison 
between many sources. 

2. OBSERVATIONS AND METHODS 

2.1. Sample and Data Analysis 

We analyzed archival Chandra ACIS observations of 
the seventeen SNRs listed in Table 1. Seven of our 
sources are thought to be Type la SNe, while ten are 
considered to have originated from a CC SNe (see Ta- 
ble 1). We selected only sources in the Milky Way 
and LMC that had their entire X-ray emitting area ob- 
served with one ACIS pointing; this criterion excluded 
sources with large spatial extent, such as SN 1006. We 
also required the targets to have strong line emission at- 
tributed to the ejecta, with Si XIII line emission surface 
brightnesses >0.01 counts/pixel 2 within the apertures R 
listed in Table 1. This criterion excluded SNRs with 
weak lines and those dominated by continuum emission 
(bremsstrahlung or non-thermal), e.g. Gl.9+0.3. This 
requirement removed two SNRs whose X-ray line emis- 
sion is attributed to shocked int erstellar medium ( ISM) 
and not e jecta: SNR 0506- 68.0 (jHughes et alf 2QQ6) and 
3C 391 (jChen et al.l 12004( 1. Additionally, we excluded 
sources with large-scale morphologies that are the re- 
sult of axisymmetric winds from pulsars (like Crab or 
3C 58). One other source, N63A, was removed since its 
morphology is distorted by its association with an HII 
region and obscuration by an intervening dense cloud 
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TABLE 1 
Archival SNRs, Sorted by Age 



Source ObsID ACIS Exp. Age Distance Radius a Lx h References 

(ks) (years) (kpc) (pc) (xlO 37 erg/s) 



Type la Sources 

0509-67.5 776, 7635, 8554 113 350-450 50 5^96 L76 Bad enes et al. (2008) 

Kepler 6714-6718, 7366 751 405 5.0 3.88 0.19 Reynolds et al. (2007) 

Tycho 3887 150 437 2.4 3.72 0.12 Warren et al. (2005) 

0519-69.0 118 40 400-800 50 6.56 1.06 Rest et al.72005) 

N103B 125 37 860 50 5.96 1.67 Lewis et al. (2003) 

DEM L71 775, 3876, 4440 148 -4360 50 11.9 0.77 Hughes et al. (2003) 

0548-70.4 1992 60 ~7100 50 17.9 0.96 Hendrick et al.(2003) 



Core-collapse Sources 



Cas A 


4634-4639, 5196, 5319-5320 


993 


309-347 


3.4 


3.77 


2.58 


Hwang et al. (2004) 


W49B 


117 


55 


~1000 


8.0 


6.30 


4.48 


Lopez et al. (2009) 


G15.9+0.2 


5530, 6288, 6289 


30 


-1000 


8.5 


8.72 


1.38 


Reynolds et al (2006) 


Gil. 2-0.3 


780-781, 2322, 3909-3912 


95 


1623 


5.0 


4.17 


1.18 


Kaspi et al. (2001) 


Kes 73 


729 


30 


500-2200 


8.0 


6.11 


0.94 


Gotthelf fc Vasisht (1997) 


RCW 103 


970 


49 


~2000 


3.3 


5.43 


2.21 


Carter et al. (1997) 


N132D 


5532, 7259, 7266 


90 


-3150 


50 


21.47 


9.92 


Borkowski et al. (2007) 


G292.0+1.8 


6677-6680, 8221, 8447 


516 


-3300 


6.0 


9.45 


0.56 


Park et al. (2007T" 


N49B 


1041 


35 


10000 


50 


23.85 


2.39 


Park et al. (2003) 


B0453-685 


1990 


40 


13000 


50 


20.27 


0.54 


Gaensler et al. (2003) 



a Radius R used in Eq. 3 of the power-ratio calculation. This radius is selected to enclose the entire source in the full-band 
X-ray image, and it is determined assuming the distances listed above. 

b X-ray luminosity in the 0.5-2.1 keV band, from the Chandra SNR Catal og. The value for G15.9+0.2 is calculated using 
an absorbed planar shock (vpshock) model of the integrated spectra (Reynolds et al. 2006). 



(IChu fe KennicuttiH98a TWarren et al.ll2QQ3h . The final 
sample reflects all SNRs that meet the above criteria. 

Each source was observed for ^30-1000 ks. Data re- 
duction and analysis was performed using the Chandra 
Interactive Analysis of Observations (ciao) Version 4.0. 
We followed the CIAO data preparation thread to repro- 
cess the Level 2 X-ray data, and we extracted the global 
Chandra X-ray spectra of each source using the CIAO 
command specextract. Exposure-corrected images of the 
Si xm line were then produced for each source by fil- 
tering data to the narrow energy range of that feature 
(~1. 75-2.00 keV). The Si xm line was selected for our 
comparative analysis as it is the only strong emission line 
common to all seventeen sources. For sources with bright 
pulsars (Gll.2-0.3, Kes 73, RCW 103, G292.0+1.8), the 
pulsar location and extent was identified using the CIAO 
command wav detect (a source detection algorithm using 
wavelet analysis; Freeman et al. 2002). We replaced 
the region identified by wavdetect with pixel count val- 
ues selected from the Poisson distribution of the area 
surrounding the pulsar using the CIAO command dmfilth. 
This process removed the bright pulsar while preserving 
the morphology of the emission surrounding the pulsar. 
No other point sources were removed because of potential 
confusion with small ejecta clumps. 

Given the young-to-middle age of our sources (see Ta- 
ble 1), we expect that all are ejecta-dominated. X-ray 
spectral modeling generally confirms that abundances 
are above those of the ISM (see references in Table 
1), suggesting the emission is indeed dominated by the 
shocked SN eje cta and not by sho cked ISM. In the case 
of G292.0+1.8. lPark et al.l (|2004T ) show that the Si xm 
and S XV is contaminated by shocked ISM in its "belt" . 
We note that our analysis of Ne IX and Mg XI in this 
source (which Park et al. 2004 attribute to ejecta) gave 



results consistent those of the Si xm line. 

2.2. Methods 

To measure the asymmetry of the X-ray line emission, 
we have applied a power-ratio method (PRM) to the Si 
xm images of our seventeen sources. The PRM was first 
applied to quantify the X-ray morphology of galaxy clus- 
ters observed with ROSAT (Buote fc T sai, 1995, 1996) 
and w ith Chandra (jJeltema et al.ll2005l ). In lLopez et al l 
(2009), we employed the method to compare the symme- 
try and elongation of line emission in the complex SNR 
W49B. We refer the reader to these papers for a detailed 
formalism; here, we give a brief overview of the technique. 

The PRM measures asymmetries in an image via cal- 
culation of the multipole moments of the X-ray surface 
brightness in a circular aperture. It is derived similarly 
to the multipole expansion of the two-dimensional grav- 
itational potential within an enclosed radius R: 



0) = -2Ga In J - 2G 

oo ^ 

X ^ mR m ( ttm COS m ^ + bm Sin ' ( 1 ) 

m=l 

where the moments a rn and b rn are 

a m (R)= [ E(f) (#O m cosra0'dV, 

JR'<R 

b m (R)= [ Z(x')(R') m smm0'd 2 x', 

JR'<R 

x' = (R f , ft), and E is the surface mass density. For our 
imaging analyses, the X-ray surface brightness replaces 
surface mass density in the power ratio calculation. 
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Fig. 1. — Left: P2/P0 versus P3/P0 for sources classified as core-collapse SNRs (blue) and the Type la SNRs (red). SNR 0548-70.4 is 
plotted in purple due to its anomalous eject a properties (see Section 3). Right: Si xm images of seven sources, plotted in the location of 
their power ratio values. Cyan Xs mark their full-band centroids, and purple dots are placed at their locations on the power-ratio plot. 
The images of RCW 103 and G292.0+1.8 have their pulsars removed using the method outlined in §2. 



The powers of the multipole expansion are obtained 
by integrating the magnitude of ^ m (the rath term in 
the multipole expansion of the potential) over a circle of 
radius R, 

Pm{R) = ^ / *m(i«, 0)*m(A, (2) 

^ Jo 

Ignoring the factor of 2G, this equation reduces to 



using the program AdaptiveBin ([Sanders fc Fab ian 2001) 
such that all zero pixels are removed to smooth out noise. 
Then, noise is added back in by taking each pixel inten- 
sity as the mean of a Poisson distribution and selecting 
randomly a new intensity from that distribution. This 
process was repeated 100 times for each Si xm image, 
creating 100 mock images per source. The 1-cr confi- 
dence limits represent the sixteenth highest and lowest 
power ratio obtained from the 100 mock images of each 



P = [a In (R)} 2 

P rn = ^J^K + bli) (3) 

The moments a m and b m (and consequently, the pow- 
ers P m ) are sensitive to the morphology of the X-ray 
surface brightness distribution, and higher-order terms 
measure asymmetries at successively smaller scales rela- 
tive to the position of the aperture center (the origin). 
To normalize with respect to flux, we divide the powers 
by Po to form the power ratios, P m /Po. Pi approaches 
zero when the origin is placed at the surface-brightness 
centroid of an image, so we have set the aperture cen- 
ter in all analyses to the full-band (0.5-8.0 keV) cen- 
troid of each remnant. In this case, morphological in- 
formation is given by the higher-order terms. P2/P0 is 
the quadrupole ratio; examples of sources that give high 
P2/P0 are those with elliptical morphologies and those 
with off-center centroids because one side is substantially 
brighter than the other. P3/P0 is the octupole ratio; ex- 
amples of sources that give high P3/P0 are those with 
deviations from mirror symmetry relative to their cen- 
troid. 

A M onte Carlo approach described in lLopez et al.l 
(2009) is used to estimate the uncertainty in the power 
ratios. Specifically, the exposure-corrected images (nor- 
malized to have units of counts) are adaptively-binned 



3. DISCUSSION 

The calculated power ratios P2 / Po and P3 / Po for the Si 
xm images of our sources are plotted in Figured] (left), 
with example Si xm images of seven sources shown in 
Figure [T] (right). The SNRs with elongated or barrel-like 
features (like W49B, Cas A, and G292.0+1.8) have the 
largest P2/P0 values, and sources that are more com- 
pact and circular (e.g., 0519—69.0 and DEM L71) have 
the lowest P2/P0. Some sources that appear spherical 
by eye (like 0509-67.5 or 0548-70.4) have large P 2 /P 
because they are brighter on one side and their full-band 
centroids have off-center positions. SNRs with obvious 
mirror asymmetries (such as N132D and N103B) produce 
large P3/P0, while sources with line emission distributed 
evenly around their full-band centroids (e.g., DEM L71 
and 0509-67.5) have small P3/P0. 

Overall, we find that the power ratios of the Type la 
SNRs are significantly different than those of the CC 
SNRs. The CC SNe have roughly an order of magnitude 
larger quadrupole power ratio: the mean P2/P0 of the 
Type la SNe is (7.7±0.7) x 10~ 7 with a standard devia- 
tion of 6.9xl0" 7 (excluding SNR 0548-70.4, see discus- 
sion below) and of the CC SNe is (80.3±2.6)xl0- 7 with 
a standard deviation of 56.1 xlO -7 . This result implies 
that the CC SNR line emission is comparatively much 
more elliptical or have more off-center centroids than 
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Type la SNRs. We find that the mean P3/P0 are also dif- 
ferent: the mean of the Type la SNe is (1.7±0.1)xl0 -7 
(excluding SNR 0548—70.4) with a standard deviation 
of 2.1 xlO- 7 and of the CC SNe is (3.5±0.3) x 1(T 7 with 
a standard deviation of 1.8xl0 -7 . This result indicates 
that the X-ray line emission in CC SNe is more mirror 
asymmetric than Type la SNe. Additionally, we find that 
the Type la SNRs with greater ellipticity or off-center 
centroids (larger P2/P0) tend to have larger mirror sym- 
metry (smaller P3/P0), and those with less mirror sym- 
metry are significantly more circular or balanced surface 
brightness (smaller P2/P0). Overall, the Type la and CC 
SNe seem to naturally form two distinct populations in 
the P2/P0 - P3/P0 plane. 

Only one source of the seventeen is an outlier in Fig- 
ured! SNR 0548-70.4. This source has two bright limbs 
as well as a large, luminous area that is off-center and 
possi bly seen in pro j ection (see the right panel of Figure 
|T]). iHendrick et al.l (|2003l ) categorized this source as a 
Type la based on its low O/Fe ratio of ~16, an inter- 
mediate value b etween the O/Fe rati o of ~0.75 expected 
for a Type la (jlwamoto et al.l fl999l ) and the O/Fe ra- 
tio of 70 predicted for a 20 M CC SNe; these authors 
present several possible explanations for an elevated oxy- 
gen abundance in a Type la SN. In our opinion, further 
study of SNR 0548—70.4 is warranted to explore its ex- 
plosion type and anomalous ejecta properties. 

By including targets from the LMC and the Milky Way, 
we aimed to limit any biases in our sample. Our SNe 
span a large range of ages and X-ray luminosities (see 
Table 1). On average, our Type la sources have younger 
estimated ages than the CC SNe, but we find no trends 
in the power ratios with age or with radius (R from Table 
1). The 0.5-2.1 keV X-ray luminosities of the Type la 
and CC SNe are approximately similar, so our analysis 
is not biased by the brightness of the sources. We note 
here that we have performed the same analyses for other 



X-ray emission lines (e.g., Ne IX, Mg XI, S xv), and all 
produce the the two distinct populations in the P2/P0 ~ 
P3/P0 plane. 

Our results show that as a class, CC SNRs are much 
more asymmetric or elliptical than Type la SNRs. These 
morphological differences reflect the distinct explosions 
and CSM structures of the progenitors of Type la and 
CC SNe. This finding is consistent with the emerging 
picture from spectropolarimetry studies that CC SNe are 
intrinsically aspherical (see Wang & Wheeler 2008 for a 
review), whereas Type la SNe are largely spherical (with 
some variation; see Kasen et al. 2009). Further obser- 
vations and analyses are necessary to elucidate whether 
the asymmetries in CC SNRs are properties determined 
by the central engine/ignition process. 

As one possible end of the stellar life cycle, SNe pro- 
vide crucial tests of stellar evolution models. However, 
in order to compare observation with theory, we need 
viable progenitor diagnostics for each SN type. During 
and immediately after explosions, many aspects of the 
physical mechanism can be explored: e.g., neutrinos, po- 
larization, radioactive decays, relic compact objects. As 
presented here, the distribution of line emission in young 
SNRs gives a complementary and unique means to probe 
the explosion mechanism many years after a SN event. 
This tantalizing possibility has always been the magnet 
of young SNR studies, and the superb spatial and spec- 
tral resolution of Chandra has enabled detailed analyses 
of the metal-rich ejecta in SNRs. The clean, quantitative 
method we introduce here is based entirely on observ- 
ables, and it represents the first standardized approach 
to type young SNRs. 

We thank the referee, Patrick Slane, for his helpful 
feedback. This work is supported by DOE SciDAC DE- 
FC02-01ER41176 (LAL and ER-R) and a National Sci- 
ence Foundation Graduate Research Fellowship (LAL). 
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